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We study the dynamics of a simple model for quantum decay, where a single state is coupled to 
a set of discrete states, the pseudo continuum, each coupled to a real continuum of states. We find 
that for constant matrix elements between the single state and the pseudo continuum the decay 
occurs via one state in a certain region of the parameters, involving the Dicke and quantum Zeno 
effects. When the matrix elements are random several cases are identified. For a pseudo continuum 
with small bandwidth there are weakly damped oscillations in the probability to be in the initial 
single state. For intermediate bandwidth one finds mesoscopic fluctuations in the probability with 
amplitude inversely proportional to the square root of the volume of the pseudo continuum space. 
They last for a long time compared to the non-random case. 

PACS numbers: 03.65.-w, 85.35.Bc, 72.15.Lh, 33.25. -|-k 



The problem of the decay of an excitation into a contin- 
uum is a fundamental problem in quantum mechanics [l[ , 
appearing in numerous fields of physics. A natural hier- 
archy of couplings occurs in many physical systems, 

For example, a spin of a nucleus may be coupled to the 
electromagnetic modes of a cavity in which it is situated, 
and these in turn may be coupled to the modes of a larger 
Hierarchical systems were studied 
but the time dependence for the 



box or the vacuum 



model we present was not investigated. 
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FIG. 1: The model consists of a single state (SS) coupled to 
a discrete number of states which form the pseudo continuum 
(PC), each of which is coupled to a dense and broad 'real' 
continuum (RC). 

Definition of the problem and main results - Our model 
is as follows: The single state (SS) |0) is coupled with ma- 
trix elements to N states which form the pseudo con- 
tinuum (PC) , of bandwidth D and constant level spacing 
d = j^zii ■ Each of these is coupled with a matrix element 

to a broad and dense real continuum (RC), see Fig. [TJ 
Notice that for each state of the PC the matrix element is 
identical for all states in the continuum, see footnote 11 1. 
For convenience we take Eq = 0, and all matrix elements 



ability amplitude of finding the system in state |0) is 
denoted by Ao{t). Starting with ^o(O) = 1, we evaluate 
Ao(i) at i > 0. The Hamiltonian of the system is: 

JV M 

fftotal = J2 ^mIm)(/"I + J2 ^k\k){k\ + ^couplings, (1) 
fj.=l k=l 

with 3> 1, and M — > oo (keeping the level spacing 
infinitesimal and the bandwidth large) and 
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We will first show that in the case of constant cou- 
plings ~ u, ~ V the system decay is dominantly 
exponential with rate: 



(3) 
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to be real, assuming time reversal symmetry. The prob- 



where 7 = -kv^Vc, being the RC density of states [12j . 
As the coupling to the RC becomes stronger, the rate de- 
creases, a phenomenon referred to as the quantum Zcno 
effect [ISl, |lj] . A diagonalization of Hamlitonian ^ for 
the case 7 ^ d yields a very wide state, similar to the 
Dicke state [l5|, and another wide state, which is the 
most relevant for the decay of the system. The other 
A^ — 1 states acquire a small width and arc not pertinent 
for the dynamics of the system until long times of order 

When the couplings to the RC are random we iden- 
tify four regimes. The boundaries between them are con- 
trolled by the typical matrix element u ^ . For < d 
there are decaying Rabi oscillations via one state. For 
uq < D Fermi's golden rule (FGR) is obtained, while 
for Uq ^ D novel mesoscopic fluctuations with amplitude 
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appear, see Fig. [S] For the almost degenerate 
case, i.e., uq > D, we find decaying oscillations via a 
linear combination of many states in the PC. 

Physical realizations.- The model ([T|) applies, for ex- 
ample, to a small quantum dot coupled to one or more 
quantum dots, each of which is coupled to a lead. One 
can study the time dependence of an injected electron's 
probability to remain in the dot, assuming that all rel- 
evant levels in the dots are empty. In some cases the 
matrix elements u^^ connecting two dots can be taken as 
constant [l^. For the more generic case of a single dot 
coupled to a disordered or sufficiently distorted larger 
quantum dot, do not have the same sign, and we take 
them to be random. These solid state implementations 
have close analogies when one replaces the quantum dots 
by atoms in optical cavities 

Derivation of the results.- Eliminating the amplitudes 
in the RC after Laplace transforming the equations of 
motion, the dynamics of the other N -\- 1 amplitudes in 
the SS and the PC are formulated in terms of a iV + 1 by 
N-\-l non-hermitian matrix, which is the matrix describ- 
ing the original matrix elements between these states, 
plus a matrix element —i^f^, = —iirVcV^i 
reduced Hamiltonian for the system is: 



(4) 

From now on we shall assume to be constant, and 
therefore 7^^ = 7. 

This leads to the (exact) eigenvalue equation: 
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The generic form for the eigenvectors is: 
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where A„ is the corresponding eigenvalue. 

Although H is non-hermitian, we can still decom pos e 
the initial state as a superposition of its eigenvectors [l8| : 



Aoit)=Y,Cn 
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Notice that even when we normalize the states the 
coefficients Cn are not the usual projections (l^n|0) [H, 

We shall now use this formalism to study the cases of 
constant and random matrix elements. 



A. Constant matrix elements.- First we take the matrix 
element between the initial SS and the levels of the PC 
to be a constant u 

To analyze the decay, let us write the equations of mo- 
tion for the amplitudes following from Eq. (|4|): 



' — 
dt 



lAq - 17 
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where is the amplitude of state p.. 

Upon Laplace transforming Eq. ([8]), we obtain: 



iw^o = i + u^A,., (9) 

iu:A^, = uAq - ^7^^^ + E^,A^,. (10) 

For D ^ uj ^ d and for the SS with energy far enough 
from the edges of the PC band, we can approximate the 



iuj — E. 



by — which leads to the result: 



An = 



r' 



(11) 



with F given by Eq. ([3]). The inverse Laplace trans- 
form gives the exponential decay, in a large time window, 
which for 7 ^ d is given by < t < ^. 

We shall now analyze the structure of the eigenstates 
yielding the result of Eq. ([3|) . 

For the limit 7 = 0, we have the discrete Wigner- 
Weisskopf problem [H, and Eq. ([3]) reduces to FGR. 
In that case the eigenvalues are real (since the Hamilto- 
nian is hermitian), and it is the superposition of many 
eigenvectors that gives rise to the decay (for intermediate 
times). 

As 7 increases and reaches the regime 7 S> d, the be- 
havior is changed and one state completely dominates the 
decay of the system. In this regime there is a fast decay- 
ing eigenvector approximately of the form |0) -f- x 

with eigenvalue A « —i^N and x « This is related 
to the Dicke effect, where a coherent sum of many states 

with equal amplitudes is also present [l^. We also find 

2 

an additional eigenvalue — i-^- Since an increase in 7 
causes a smaller decay rate, we are motivated to call the 
corresponding eigenvector the Zeno state. Using pertur- 
bation theory [18| one may show that in this case the 
coefficients C„ of Eq. ([7]) are nearly unity for the Zeno 
state, and much smaller than unity for all other states. 
Since the other states decay much slower, at long times 
the Zeno state stops being the dominant state in the de- 
cay process. Their weight in the decomposition can be 
bounded by ^ , and using the Zeno state decay rate yields 
a crossover time of order :^log(^) [IS]- The perturba- 
tion theory results are confirmed numerically, see Fig. [21 
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FIG. 2: Imaginary part of the eigenvalues of Hamiltonian (j4|, 
with constant matrix elements, u is 0.1, D is 1, 7 is 1, and 
there are 101 states in the PC. The inset shows \C„ \ define in 
Eq. (O. The weight is concentrated on one eigenvector, the 
Zeno state, decaying faster than all the rest of the eigenvalues 
except for the Dicke state. 



B. Random matrix elements - When the disorder is suf- 
ficiently large, the matrix elements can be considered ran- 
dom |10| . For simplicity, let us consider the case where 
the level spacing is constant, but the elements are ran- 
domly distributed around 0, with a standard deviation u. 

To understand the magnitudes of the matrix elements 
involved in the physical realization of a small quantum 
dot coupled to a larger disordered one which is coupled 
to a one channel lead , it is instructive to look at their 
site representation. If we denote the sites of the larger 
dot by \i), then the isolated dot eigenstates are = 



where 0^ are random coefficients of order 
unity (assuming the disorder is large enough, yet not too 
large as to make the states localized). We shall assume 
that out of the sites, S are coupled to |0), and S' are 
coupled to the RC. Changing basis to the set one can 
verify that the couplings to the RC and |0) are random. If 
the matrix elements Vj^k (between a site j and a state k in 
the lead) do not depend on j, after the lead elimination 
the Hamiltonian contains terms ~ \^){v\af^a^^ and by 
multiplying the states by a phase factor we can obtain 
a model corresponding to Eq. ([4]), with 7^1^ real and 
positive. For simplicity we shall assume the magnitude 
of 7^1, to also be a constant, 7 l2l|. Denoting the typical 
matrix element connecting a site in the larger dot with 
a site in the lead by vq and the typical matrix element 
connecting a site in the larger dot with a site in the single 
state dot uq , a straightforward calculation shows that the 
typical tunneling matrix elements in the Hamiltonian (j4|) 
are given by the relations: 7 = 7o-^: u ~ "^'"o where 
70 = ^'cl^^oP and u is the typical matrix element in Eq. 
([j]). Notice that uq, fOj7o are microscopic parameters 



independent of the system size. 

We shall now analyze the dynamics for four different 
cases: uq > D , uq ^ D , D > uq > d, d > uq. 

Case I.- To understand the behavior for uq > D, we 
study the degenerate case D = 0, for which the Hamilto- 
nian is H = E^«,.(|0)(mI + Im)(0|) - *7Em,. \f^)H 

Defining the states \Wi) = ^==E and 

IVF3) = 1/^): the Hamiltonian takes the form: 



H = 



+ - ^N^\W3){W3\. (12) 



As iM^i) and jM^a) are not orthogonal, it is useful 
to use a Gram-Schmidt procedure to define 111^2) = 
\Wi) — |II^3)(Wi|W3) and then normalize it. Now we can 
represent the system's Hamiltonian in the basis formed 
by |0), 114^2) and {Ws) as a 3x3 matrix: 



U2^{l-C^) CU2 

H= \ U^Vi^-c') 

CU2 -ijN 



(13) 



where Ui 

Ui 1 



U2VN Vn' 

Diagonalizing H perturbativley in we find (using 
c ~ --^) eigenvectors IV3) = IW3) + 0(-^) with eigen- 
value -i-/o + 0{j^),andV± = ^(|0) ± jl^a)) with eigen- 
values A± = ±1/2 ~ iS, where 5 = -i^^ ^T^- 

Since the projection of the initial state on | V3) is neg- 
ligible, the decay of the initial state is described by a 
superposition of two exponentially decaying terms, with 
exponents « ±J72 — i5. This implies Rabi-type oscilla- 
tions ([H, vol 1, p. 447) with a characteristic frequency 
U2 ^ itQ, and with an envelope decaying exponentially 



with rate 



Notice that in the case of fi- 
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nite bandwidth, the above analysis will be approximately 
correct if | A± | S> £> and 7 ^ d . The first condition gives 
the restriction D <C uq- Since this relation is N indepen- 
dent, both cases are physically accessible in the limit of 
large N . The result of a numerical simulation is shown 
in Fig. Era). 

Notice two peculiarities of the result: First, as in 
Eq. ([3]), when 7 is increased the decay is slower. Second, 
the decay is much slower than in the ordered case, by a 
factor of 2N . This can be understood as follows: the sum 
of all imaginary parts of the eigenvalues exactly equals 
—iNj. Since the Dicke eigenvector is still of the same 
form as before, with an eigenvalue approximately given 
by —iN'y, the rest of the eigenvalues have small imagi- 
nary parts. Plugging the Dicke eigenvalue into Eq. ([5]) as 

—iNj + e, one finds that e « 7^"^ . For constant ma- 
trix elements, the part of the decay sum rule not taken up 
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FIG. 3: Decay of a SS coupled to a PC with random matrix 
elements (see Eq. Q), uniformly distributed with stan- 
dard deviation 0.1. For L» = (left) 7 = 0.1, = 21, while 
for the non degenerate case D = 1 (right) 7 = 1, = 101. 
The degenerate case shows damped oscillations, while the fi- 
nite bandwidth case shows a sharp initial decay followed by 
mesoscopic fluctuations. The inset shows the short-time fit to 
a sine function. 

2 

by the Dicke eigenvector, namely was mostly taken 
up by the Zeno eigenvector. Here, the sum of all other 

-2 

decay rates ~ distributed among all eigenvalues. 

Case II.- For D ^ uq the 'disorder' is not 'self- 
averaging', and the decay depends on the realization of 
the disorder. An example is shown in Fig. [3jb). The 
initial smooth decay is approximately a sine (which, in 
the general case, will be replaced by the Fourier trans- 
form of the density of states), as shown in the inset, up to 
times of order in which the amplitude decays almost 
to 0. Then random oscillations come into play, with an 
amplitude of order . 

One can verify that the real parts of the eigenval- 
ues form a band of bandwidth D with constant density, 
by considering the eigenvalue Eq. ([5]). If Si happens 
to be exactly imaginary, then it is easy to see graphi- 
cally that there is a solution A between each consecutive 
energies i?^. Using Eq. ([6]), one can now check self- 
consistently, that Si is almost purely imaginary, since 
7 S \-E ^ ■ Therefore the eigenvalues are almost 
real, and are distributed with constant as claimed above. 
Since the eigenvectors are nearly real, so are the coeffi- 
cients in Eq. JT]), but these are now randomly distributed. 
We have the interference of N oscillating components 
with positive, random coefficients. Then the amplitude 
is described by Ao{t) « g^e"*""*, where qn arc inde- 
pendent, random variables. To understand the behavior 
in this case, it is instructive to look at the ensemble av- 
erage of probability, denoted by an overlinc. This gives: 

Performing the sum under the assumption id <C 1 gives 
|Ao(t)P ~ sinc(Dt), which is demonstrated in Fig. [3jb). 
Although we have random oscillations, a characteristic 



frequency seems apparent. This is because typically a 
few states are, by chance, coupled more than the others, 
bmce m our case q„ ^ |(y„|0)| ^ the amplitude of 

the oscillations ~ — ^ as well. 

VN 

Case III- Defining Aq = ttu^ /d ^ u^/D, the FGR rate, 
we find different behavior for the cases D ^ Aq. Notice 
that the condition of the crossover D ^ Ao is the same 
as before, namely D ^ uq. 

li D ^ Aq {D S> uq), a coefficient of an eigenvalue 
A will have a typical size Y+^TZTJa- Notice that this fits 

the exact sum ^ = 1. Thus, we have a superposition 
of M ~ oscillatory signals, with (positive) random 
coefficients. If furthermore u ^ d (which is equivalent 
to uq^^ S> D), we will have a large number of random 
components within a complete Lorentzian, and therefore 
the FGR exponential decay will be retrieved (but with 
the Zeno effect suppressed). 

Case IV.- For ito <C d, a single state is relevant, and 
the system will show (decaying) Rabi oscillations. 

Conclusions.- We considered the generic problem of a 
single state coupled to a real continuum via a pseudo 
continuum. In the ordered case, we found that a single 
eigenvector characterizes most of the decay of the sys- 
tem, and the decay becomes slower when increasing the 
coupling to the real continuum. When the bandwidth 
D is smaller than the typical matrix element uq, adding 
disorder causes the decay to be much slower, and intro- 
duces oscillations in time. When D ^ uq, mesoscopic 
fluctuations in the probability to stay in the single state 
as a function of time follow, while for a larger bandwidth, 
FGR exponential decay is retrieved. 
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